Search for: All records

Abstract The advancements in high-throughput technologies provide exciting opportunities to obtain multi-omics data from the same individuals in a biomedical study, and joint analyses of data from multiple sources offer many benefits. However, the occurrence of missing values is an inevitable issue in multi-omics data because measurements such as mRNA gene expression levels often require invasive tissue sampling from patients. Common approaches for addressing missing measurements include analyses based on observations with complete data or multiple imputation methods. In this paper, we propose a novel integrative multi-omics analytical framework based onp-value weight adjustment in order to incorporate observations with incomplete data into the analysis. By splitting the data into a complete set with full information and an incomplete set with missing measurements, we introduce mechanisms to derive weights and weight-adjustedp-values from the two sets. Through simulation analyses, we demonstrate that the proposed framework achieves considerable statistical power gains compared to a complete case analysis or multiple imputation approaches. We illustrate the implementation of our proposed framework in a study of preterm infant birth weights by a joint analysis of DNA methylation, mRNA, and the phenotypic outcome. Supplementary materials accompanying this paper appear online. 

Abstract Mixing along isopycnals plays an important role in the transport and uptake of oceanic tracers. Isopycnal mixing is commonly quantified by a tracer diffusivity. Previous studies have estimated the tracer diffusivity using the rate of dispersion of surface drifters, subsurface floats, or numerical particles advected by satellite‐derived velocity fields. This study shows that the diffusivity can be more efficiently estimated from the dispersion of coherent mesoscale eddies. Coherent eddies are identified and tracked as the persistent sea surface height extrema in both a two‐layer quasigeostrophic (QG) model and an idealized primitive equation (PE) model. The Lagrangian diffusivity is estimated using the tracks of these coherent eddies and compared to the diagnosed Eulerian diffusivity. It is found that the meridional coherent eddy diffusivity approaches a stable value within about 20–40 days in both models. In the QG model, the coherent eddy diffusivity is a good approximation to the upper‐layer tracer diffusivity in a broad range of flow regimes, except for small values of bottom friction or planetary vorticity gradient, where the motions of same‐sign eddies are correlated over long distances. In the PE model, the tracer diffusivity has a complicated vertical structure and the coherent eddy diffusivity is correlated with the tracer diffusivity at the e‐folding depth of the energy‐containing eddies where the intrinsic speed of the coherent eddies matches the rms eddy velocity. These results suggest that the oceanic tracer diffusivity at depth can be estimated from the movements of coherent mesoscale eddies, which are routinely tracked from satellite observations. 

Abstract The vertical structure of ocean eddies is generally surface-intensified, commonly attributed to the dominant baroclinic modes arising from the boundary conditions (BCs). Conventional BC considerations mostly focus on either flat- or rough-bottom conditions. The impact of surface buoyancy anomalies—often represented by surface potential vorticity (PV) anomalies—has not been fully explored. Here, we study the role of the surface PV in setting the vertical distribution of eddy kinetic energy (EKE) in an idealized adiabatic ocean model driven by wind stress. The simulated EKE profile in the extratropical ocean tends to peak at the surface and have ane-folding depth typically smaller than half of the ocean depth. This vertical structure can be reasonably represented by a single surface quasigeostrophic (SQG) mode at the energy-containing scale resulting from the large-scale PV structure. Due to isopycnal outcropping and interior PV homogenization, the surface meridional PV gradient is substantially stronger than the interior PV gradient, yielding surface-trapped baroclinically unstable modes with horizontal scales comparable to or smaller than the deformation radius. These surface-trapped eddies then grow in size both horizontally and vertically through an inverse energy cascade up to the energy-containing scale, which dominates the vertical distribution of EKE. As for smaller horizontal scales, the EKE distribution decays faster with depth. Guided by this interpretation, an SQG-based scale-aware parameterization of the EKE profile is proposed. Preliminary offline diagnosis of a high-resolution simulation shows the proposed scheme successfully reproducing the dependence of the vertical structure of EKE on the horizontal grid resolution. 

Abstract Isopycnal mixing of tracers is important for ocean dynamics and biogeochemistry. Previous studies have primarily focused on the horizontal structure of mixing, but what controls its vertical structure is still unclear. This study investigates the vertical structure of the isopycnal tracer diffusivity diagnosed by a multiple‐tracer inversion method in an idealized basin circulation model. The first two eigenvalues of the symmetric part of the 3D diffusivity tensor are approximately tangent to isopycnal surfaces. The isopycnal mixing is anisotropic, with principal directions of the large and small diffusivities generally oriented along and across the mean flow direction. The cross‐stream diffusivity can be reconstructed from the along‐stream diffusivity after accounting for suppression of mixing by the mean flow. In the circumpolar channel and the upper ocean in the gyres, the vertical structure of the along‐stream diffusivity follows that of the rms eddy velocity times a depth‐independent local energy‐containing scale estimated from the sea surface height. The diffusivity in the deep ocean in the gyres instead follows the profile of the eddy kinetic energy times a depth‐independent mixing time scale. The transition between the two mixing regimes is attributed to the dominance of nonlinear interactions and linear waves in the upper and deep ocean, respectively, distinguished by a nonlinearity parameter. A formula is proposed that accounts for both regimes and captures the vertical variation of diffusivities better than extant theories. These results inform efforts to parameterize the vertical structure of isopycnal mixing in coarse‐resolution ocean models. 

x-ray polarization measurements determine the geometric arrangement of hot material accreting onto a black hole.